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ABSTRACT 

We present mid-infrared spectroscopy obtained with the Spitzer Space Telescope of a sample of 1 1 optically 
faint, infrared luminous galaxies selected from a Spitzer MIPS 70/im imaging survey of the NDWFS Bootes 
field. These are the first Spitzer IRS spectra presented of distant 70/j,m-selected sources. All the galaxies lie 
at redshifts 0.3 < z < L3 and have very large infrared luminosities of L/r ^^O.l-Hx 10'^ Lq. Seven of the 
galaxies exhibit strong emission features attributed to polycyclic aromatic hydrocarbons (PAHs). The average 
IRS spectrum of these sources is characteristic of classical starburst galaxies, but with much larger infrared 
luminosities. The PAH luminosities of vh^il .1 pva)^ 0.4-7 x 10" L© imply star formation rates of ~ 40-720 
Mq yr"'. Four of the galaxies show deep 9.1 pm silicate absorption features and no significant PAH emission 
features (6.2/xm equivalent widths < 0.03/xm). The large infrared luminosities and low vfy(lQ^m)lvfy{2Apm) 
flux density ratios suggests that these sources have AGN as the dominant origin of their large mid-infrared 
luminosities, although deeply embedded but luminous starbursts cannot be ruled out. If the absorbed sources 
are AGN-dominated, a significant fraction of all far-infrared bright, optically faint sources may be dominated 
by AGN. 

Subject headings: galaxies: active — galaxies: starburst — infrared: galaxies — quasars: general 



1. INTRODUCTION 

Luminous and ultra-luminous infrared galaxies (LIRGs: 
L(8-1000/im)> 10" Lq; Sanders etal. 1988) have been stud- 
ied extensively in the local Universe with the Infrared Astro- 
nomical Telescope (IRAS; Soifer et al. 1987; Saunders et al. 
1990), the Infrared Space Observatory (ISO; e.g., Lutz et al. 
1998; Genzel & Cesarsky 2000; Tran et al. 2001), and more 
recently, with the Infrared Spectrograph (IRS; Houck et al. 
2004) on Spitzer (Weedman et al. 2005; Brandl et al. 2006; 
Armus et al. 2007; Desai et al. 2007). These galaxies ex- 
hibit a large range of properties in the mid-IR, some showing 
strong PAH emission features characteristic of powerful (up 
to RslOOO Mq yr~') star formation rates (e.g., Brandl et al. 
2006; Smith et al. 2007), and all exhibiting a large range in 
9.7/im silicate absorption or emission strengths (e.g.. Weed- 
man et al. 2005; Desai et al. 2007; Imanishi et al. 2007). 
Spitzer IRS is now enabling the study of the mid-infrared 
spectra of LIRGs to much higher redshifts (z ~2.6; Houck 
et al. 2005 Yan et al. 2005). Although rare locally, LIRGs be- 
come an important population at high redshifts and account 
for an increasing fraction of the star-formation activity in the 
Universe (Le Floc'h et al. 2005). By studying their infrared 
properties, one can estimate the extent to which AGN and 
star-formation contribute to their infrared luminosities, and 
therefore determine a correct census of starbursts and AGN at 
epochs in the Universe when their luminosity density was at 
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its maximum. 

A particularly interesting population of high redshift op- 
tically faint infrared sources has been discovered using the 
Multiband Imaging Photometer (MIPS; Rieke et al. 2004). 
Various observing programs with Spitzer IRS have found that 
MIPS sources at flux density levels of fv{2Apm) ^ 1 mJy 
with optical magnitudes > 24 Vega magnitudes are typi- 
cally at z 2 (Houck et al. 2005; Yan et al. 2007; Weedman 
et al. 2006b). Most of the sources are characterized by strong 
absorption by the 9.1pm silicate feature, but sources chosen 
with an additional indicator of star formation (sub-mm detec- 
tion or shape of the spectral energy distribution) often shown 
strong polycyclic aromatic hydrocarbon (PAH) emission fea- 
tures (Lutz et al. 2005; Weedman et al. 2006a; Menendez- 
Delmestre et al. 2007; Yan et al. 2007; Sajina et al. 2007). Us- 
ing traditional optical techniques to characterize these sources 
is difficult because they are too faint in the optical band. Mid- 
infrared spectroscopy is the best currently available tool in 
understanding these sources but we need to learn how to cat- 
egorize them from the different features that are exhibited in 
the infrared band. We might expect sources which are se- 
lected on the basis of their large far-infrared flux densities to 
contain large amounts of cool dust and show signatures of 
starbursts. Mid-infrared spectroscopy of these sources will 
test this assumption and may help our understanding of high 
redshift LIRGs as a whole. 

In this paper, we present the first Spitzer IRS spectra of a 
small sample of galaxies selected on the basis of their large 
70/im flux densities and optical faintness. We have previously 
reported results on sources chosen only with criteria of large 
infrared to optical flux ratios, as measured by comparison of 
their 24 fim flux density, f,y(24/im) with optical magnitude {R- 
[24]>14; Houck et al. 2005; Weedman et al. 2006b). We now 
extend the IRS observations to additional sources in Bootes 
that are less extreme in their /?-[24] colors but are selected 
with the additional criteria of detection at 70/im by Spitzer 
MIPS. 

A cosmology of Hq = 70kms~' Mpc"^ $7^=0. 3, and 
J7a=0.7 is assumed throughout. 



2. SOURCE SELECTION 

We selected a sample of optically faint far-infrared lumi- 
nous sources from the Spitzer MIPS 70/im survey of the 
NOAO Deep Wide-Field (NDWFS; Jannuzi & Dey 1999) 
Bootes field for follow-up with Spitzer IRS. Such sources 
have not previously been studied and our paper provides the 
first census of their basic mid-IR properties. The 70/im ob- 
servations reach a 5a limiting flux density of f70/im=25 mJy, 
and yield a total of « 330 sources. To insure a highly reliable 
70/Ltm catalog, and allow follow-up observations with IRS, 
we only included 70yKm sources with a 24/xm flux density, 
f24 >1 mJy and fjo^m >30 mJy (this resulted in only 16 re- 
jected sources and should not bias our sample to sources with 
unusually high 24/im flux densities). We measured optical 
photometry in the publicly available NDWFS /?-band images, 
and selected all sources with R >20 Vega mag. Although ide- 
ally, we would like to select sources as red in their optical 
to infrared colors as the /?-[24]>14 sources previously ob- 
served in the Bootes field, the lack of available candidates has 
required us to relax this criteria. Our sample have a larger 
range of optical to infrared colors corresponding to 11 < R- 
[24]<16.5. Although the optical obscuration is likely to be 
less extreme, the combined faint optical and bright infrared 
selection should still select distant infrared bright galaxies. 
The basic optical and infrared properties of the 1 1 sources are 
presented in Table 1 . The IRAC photometry is from the IRAC 
Shallow Survey (Eisenhardt et al. 2004). 

3. OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations were made with the IRS 
Short Low module in order 1 only (SLl) and with the Long 
Low module in orders 1 and 2 (LLl and LL2), described in 
Houck et al. (2004). These orders give low resolution spectral 
coverage from ~8 /xm to ~35 fim. Sources were placed on the 
slits by offsetting from a nearby 2MASS star. The integration 
times for individual sources are given in Table 2. 

Because these faint sources are dominated by background 
signal, we restrict the number of pixels used to define the 
source spectrum, applying an average extraction width of only 
4 pixels (which scales with wavelength). This improves the 
signal-to-noise ratio, although some source flux in outlying 
pixels is lost so a correction is needed to change the fluxes 
obtained with the narrow extraction to the fluxes that would 
be measured with a standard extraction. This flux correction 
is derived empirically by extracting an unresolved source of 
high S/N with both techniques and is a correction of about 
10%, although the correction varies with order and with wave- 
length. 

The background which was subtracted for LLl or LL2 in- 
cludes co-added backgrounds from both nod positions having 
the source in the other slit (i.e., both nods on the LLl slit 
when the source is in the LL2 slit), added together with the 
alternative nod position in the same slit, yields a background 
observation with three times the integration time as for the 
source. For SLl, there was no separate background obser- 
vation with the source in the SL2 sUt, so background sub- 
traction was done between co-added images of the two nod 
positions in SLl. Independent spectral extractions for each 
nod position were compared to reject highly outlying pixels 
in either spectrum, and a final mean spectrum was produced. 
Data were processed with version 13.0 of the SSC pipeUne, 
and extraction of source spectra was done with the SMART 
analysis package (Higdon et al. 2004). 



4. IRS SPECTRAL CHARACTERISTICS 

The IRS spectra of the eleven 70/im selected sources are 
presented in Figure 1 . They are boxcar smoothed over a reso- 
lution element (approximately two pixels). Most of the sam- 
ple (7 of 11) exhibit PAH emission features with measurable 
6.2/tm PAH equivalent widths in their infrared spectra and are 
hereafter refered to as PAH-dominated sources (denoted as 
"sb" in Table 2). The large PAH equivalent widths are consis- 
tent with them being starburst-dominated sources and was an- 
ticipated because the large 70/im fluxes imply that the spectra 
are dominated by dust that is cooler than that typically asso- 
ciated with AGN. The remaining four sources exhibit strong 
9.7 /im silicate absorption features, no obvious PAH emission 
features, and are classified as absorption-dominated sources 
(denoted as "abs" in Table 2). 

For the PAH-dominated sources, redshifts are determined 
from the strong PAH emission features at rest wavelengths 
6.2/im, 7.7/im, and 11.3/im. With sufficient signal-to-noise 
ratios, use of these features gives redshifts consistent with 
optically-derived redshifts to ±0.001 (Houck et al. 2006). 
Redshift precision for faint sources depends on the signal- 
to-noise ratio, so for the spectra with PAH features, we con- 
sider the redshifts accurate only to ±0.05. For the absorbed 
sources, redshifts are determined primarily from the localized 
maximum in the IRS continuum, or "hump", which is blue- 
ward of the absorption feature. In the absorption-dominated 
sources, this feature is produced by absorption on either side 
of the feature, and with a possible contribution of 7.7/im PAH 
emission in composite sources. We determine redshifts for 
absorption-dominated sources by assuming this hump has a 
rest-frame wavelength of 7.9/im (see average spectra in Hao 
et al. (2006) and Spoon et al. (2007)). Comparison with op- 
tical redshifts (e.g.. Brand et al. 2007) suggests that they can 
be uncertain to ±0.1 in z. 

We constructed an average mid-infrared spectra of the 7 
sources classified as PAH-dominated and for the 4 sources 
classified as absorption-dominated. This was achieved by 
wavelength correcting each individual spectrum to the rest 
frame (z = 0), interpolating them to a common wavelength 
scale of wO.l/tm/pix, and taking a straight (non-weighted) 
average at each wavelength position. The resulting average 
spectra are shown in Figure 2. We discuss each class in turn. 



4.1. PAH-dominated 70 /im sources 

The average IRS spectrum of the 7 PAH-dominated sources 
in Figure 2 is characteristic of classical starbursts. The spec- 
trum is very similar in shape to the average low redshift (and 
lower luminosity) starburst template presented by Brandl et al. 
(2006), although the PAH emission features at 11.3/im and 
12.7/im appear to be weaker in relation to the 7.7/im PAH 
emission feature. The presence of low ionization [Ne llj 
A12.81/im and [Ne III] A15.56/im emission lines and lack of 
significant high ionization features such as [Ne V] A14.32/tm, 
are also very typical of starbursts. There is a possible detec- 
tion of the rotational H2 S(3) A9.7/im emission line in the av- 
eraged spectrum. This feature is seen in both high and low res- 
olution IRS spectra of nearby starburst galaxies (e.g., Brandl 
et al. 2006; Higdon et al. 2006; Farrah et al. 2007; Armus 
et al. 2007) and suggests the presence of warm molecular gas. 
We measure the rest-frame 6.2/im PAH equivalent widths by 
fitting a single Gaussian and using the continuum adjacent to 
the feature, between 5.5/im and 7/tm. The values shown in 



Spitzer IRS Spectra of 70 micron Sources 



3 



TABLE 1 
Optical and IR properties. 



IRS ID 








ja.b 




13 .6 


14.5 


15. 8 


18 


124 




tl60 






name 


map 


map 


map 


map 


mjy 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


map 


70Bootesl 


SST24 J 14265 1 .9+343 1 35 


22.7 


20.6 


19.8 


16.7« 


0.11 


0.09 


0.08 


0.17 


1.39 


32.8±4.2 


<120 


11.3 


70Bootes2 


SST24 J142732. 9+324542 


22.7 


20.7 


19.9 




0.07 


0.07 


0.09 


0.21 


1.22 


34.6±4.1 


<120 


11.2 


70Bootes3 


SST24 1143639.0+345222 


23.2 


22.0 


21.1 




0.13 


0.10 


0.10 


0.12 


1.26 


35.0±6.3 


145±29 


12.6 


70Boote.s4 


SST24 1143218.1+341300 


23.3 


22.0 


21.0 


17.8'' 


0.10 


0.09 


0.09 


0.12 


1.22 


36.6±5.1 


100±20 


12.6 


70Boote.s5 


SST24 1143050.8+344848 




24.6 


22.4 




0.04 


0.04 


0.08 


0.24 


4.25 


43.2±5.1 


70±14 


16.5 


70Bootes6 


SST24 1143830.6+344412 


23.8 


21.1 


18.7 




0.13 


0.14 


0.22 


0.78 


3.19 


45.2±4.0 


<120 


12.7 


70Bootes7 


SST24 1143151.8+324327 


22.3 


20.5 


19.7 




0.20 


0.16 


0.17 


0.37 


2.18 


51.8±4.7 


135±27 


11.7 


70Bootes8 


SST24 1143341.9+330136 


22.7 


21.1 


20.2 


17.8'' 


0.08 


0.06 


0.09 


0.21 


4.72 


63.2±3.9 


105±21 


13.1 


70Bootes9 


SST24 1143820.7+340233 


22.7 


20.3 


19.3 


15.7'' 


0.40 


0.27 


0.32 


0.34 


3.46 


67.2±3.1 


245±49 


12.0 


70BooteslO 


SST24J143449.3+341014 


23.6 


21.5 


20.7 


17.5'' 


0.25 


0.23 


0.68 


0.97 


2.19 


94.5±4.5 


120±24 


12.7 


70Bootesll^ 


SST24J143205.6+325835 


21.3 


20.1 


19.8 




0.08 


0.13 


0.27 


1.3 


16.66 


115.8±12.8 


<120 


13.5 



AH quoted magnitudes are Vega magnitudes from ttie NDWFS DR3 (Jannuzi et al. in preparation). 
'^Errors on Sn , R, /, and AT-band magnitudes are <0.1. 
'^Errors on fj.e, f4.5, fg.Si and f24 flux densities are <0.l mJy. 
■'/sTj-mag from FLAMEX survey (Elston et al. 2006). 
'AT-mag from NDWFS survey Dey et al. in preparation. 
^Object 8 in Houck et al. 2007. 



TABLE 2 

IRS PROPERTIES AND DERIVED PHYSICAL PROPERTIES. 



IRS 


exp. time" 


IRS 




i/L^(7.7/tm)* 




L/r(24+70+160)'' 


SFR"" 


f(15)/f(6) 




PAHEW« 


IRS 


ID 


s 


z 


log(erg s"') 


log(erg s-i) 


log(L0) 


log(Lo) 


M© yf 


-1 




jjm 


class* 


70Bootesl 


480,1200 


0.501 


43.89 


44.53 


11.72 


12.15 


90 


4.2 




0.26±0.05 


sb 


70Bootes2 


480,1200 


0.366 


43.65 


44.18 


11.36 


11.79 


40 


1.8 




0.85±0.17 


sb 


70Bootes3 


600,1440 


0.986 


44.69 


45.37 


12.56 


13.05 


625 


4.3 




0.48±0.05 


sb 


70Bootes4 


480,1200 


0.975 


44.24 


45.32 


12.51 


12.95 


565 


9.6 




0.59±0.07 


sb 


70Bootes5 


240,480 


1.21 


45.51 






13.24 




7.1 


-1.0 


<0.02 


abs 


70Bootes6 


240,480 


0.94 


45.27 






13.01 




4.0 


-1.9 


<0.01 


abs 


70Bootes7 


480,960 


0.664 


44.46 


45.20 


12.39 


12.61 


430 


3.4 




0.29±0.04 


sb 


70Bootes8 


240,480 


0.81 


44.79 






12.88 




7.7 


-3.6 


<0.03 


abs 


70Bootes9 


240,480 


0.668 


44.62 


45.43 


12.62 


12.81 


720 


4.0 




0.53±0.04 


sb 


70BooteslO 


480,960 


0.512 


44.43 


44.80 


11.98 


12.41 


170 


2.1 




0.13±0.01 


sb 


70Bootesll' 


240,480 


0.48 


44.75 






12.49 




8.3 


-2.2 


<0.02 


abs 



^Total integration times in SLl (first entry) and each of LLl and LL2 (second entry; same for each). 

''7.7/im luminosity in source rest frame, determined from peak flux density of 7.7/im feature without continuum subtraction. 
'Infrared luminosity derived using the relation: log[L/s]=log[i'L„(7.7/:im)] + 0.78 from Houck et al. (2007). 
'^Minimum infrared luminosity estimated from L/s=[i'L^(24(Um)+i/L,y (70/:im)+i'L^(160([im)] X (1+z). 
■^Star formation rate derived from infrared Itxminosity using relation from Kennicutt (1998). 
'Sihcate absorption sttength as defined in Section, 4,2. 
^Rest-frame eqixivalent width of (>.2nm PAH emission feature, 

''PAH-dominated and absorption-dominated soixrces are denoted "sb" and "abs" respectively. See Section 4 for definitions of these terms. 
'Object 8 in Houck et al. 2007. 



Table 2 are similar to that of the lower redshift starburst galax- 
ies presented by Brandl et al. (2006) and Desai et al. (2007). 

These sources are interesting because of their extreme mid- 
infrared luminosities and red colors. The median continuum 
luminosity is z^Ly(6/im)=4 x lO"*^ ergs s~' (1 x 10' ' L©). 
For comparison, the median luminosity for starbursts in the 
Bootes f24 >10 mJy sample of Houck et al. (2007) is 1 x 
lO^"' ergs s~' (3 x 10^ L©) and for the prototype starburst 
NGC 7714 of Brandl et al. (2006) is 5 x 10^2 ergs s"' (1 
x 10^ Lq). In Figure 3 we compare the PAH luminosities, 
\og{vL .1 jim)), of our sample to other starbursts. The mean 
PAH luminosity of the PAH-dominated sources is orders of 
magnitude larger than that of the local starbursts from Brandl 



et al. (2006). The least luminous sources in our sample are 
comparable to the most luminous starbursts in the local Uni- 
verse. The luminosities are slightly lower than that of the 
high redshift starbursts from Weedman et al. (2006a) which 
were chosen from the existence of a stellar photospheric peak 
in their near-IR SEDs. If we use the conversion in Houck 
et al. (2007) between the 7.7 /im PAH luminosity and star 
formation rate, we find that our most luminous source has a 
star-formation rate of ~720 M© yr"'. This is only ~50% 
lower than the upper Umit on the Ha estimated star forma- 
tion rates in local bright galaxies (Kennicutt 1998). Table 2 
shows that these sources have larger infrared luminosities (as 
estimated from their 24/xm, 70/xm, and 160/xm luminosities) 
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Fig. 1 . — IRS spectra of 70 um selected luminous infrared galaxies. The sources are divided into PAH-dominated (top 7) and absorption-dominated (bottom 4) 
sources. The PAH-dominated sources are ordered by their 7.7/im PAH luminosities (from strongest to weakest). The absorption-dominated sources are ordered 
by their silicate absorption strengths (shallow to deep). The spectra are boxcar smoothed over a resolution element (approximately two pixels). The expected 
positions of the PAH emission features and silicate absorption features are shown. The measured equivalent widths of (he PAH features and silicate absorpion 
depths are given in Table 2. Also plotted on the bottom right is a typical noise spectrum (for 70Bootes9). 
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Fig. 2. — Average IRS spectrum of PAH-dominated sources (top) and 
absorption-dominated sources (bottom). The expected positions of typically 
strong emission lines and the silicate absorption line are plotted. The average 
IRS spectrum of low redshift starburst galaxies from Brandl et al. (2006) is 
over-plotted for the PAH-dominated spectrum (normalized to the 7.7/jm peak 
flux; dotted line). 

than is estimated from their 7.7 /im PAH luminosities. Their 
infrared luminosities are among the largest known for star- 
forming galaxies. We measure the rest-frame fy(6/im) and 
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Fig. 3. — Histogram showing the 7.7 fim PAH luminosities, 
log(uhu(7 J /im)), of the 70 fim PAH-dominated sources in our sample (black 
filled histogram). For comparison, we show the local starbursts from Brandl 
et al. 2006 (red vertical stripes), the 10 mjy PAH-dominated sample from 
Houck et al. 2007 (blue diagonal stripes), and the high redshift starburst sam- 
ple from Weedman et al. 2006a (green horizontal stripes). 

fy(15/^m) flux densities from the IRS spectra. fy(6/im) is 
measured just shortward of 6/im to avoid the 6.2/im PAH fea- 
ture and provide a measurement of the hot dust continuum. 



In Figure 4, we show the distribution of 6/im luminosities 
(i^Le^m) as a function of continuum slope as measured by the 
rest-frame flux density ratio, i^fy(15/im)/i^f,y(6/im) and ob- 
served flux density ratio, i^fy(70/im)/j/f^(24/im). The PAH- 
dominated sources have a lower mean t^fy(15/im)/j/fy(6/im) 
than that of the lower luminosity starbursts of Brandl et al. 
(2006). Although this could suggest a contribution from an 
AGN to the infrared emission (since the hot dust results in 
a shallower spectral slope), a larger sample size is needed to 
show this with any confidence. 
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Fig. 4. — hog(uhi/(6fim)) as a function of rest- 

frame uf,^(l5fj,m)/ui,y(6fj.m) ratio (top) and observed frame 
!^f^(70/im)/i/fy(24/im) ratio (bottom). PAH-dominated sources are 
represented by large black squares; absorption-dominated sources are 
represented by large black triangles; small filled circles represent PAH 
emission sources from the Bootes 10 mJy sample (Houck et al. 2007); small 
empty circles represent local starbursts from Brandl et al. (2006); large green 
triangles represent X-ray-selected AGN from Brand et al. (2007). 



4.2. Absorption-dominated 70 fim sources 

The four absorption-dominated 70/im sources have IRS 
spectra which are generally more characteristic of AGN: they 
have no obvious PAH emission features (6.2/im PAH equiv- 
alent width < 0.03/im) but deep silicate absorption features. 
They have 6/im luminosities that are larger than that of the 
PAH-dominated sources and similar to that of the X-ray se- 
lected AGN in Brand et al. (2007). The infrared luminosities 
are also very large (L/r ^ 0.3 -1.7x10'-' L©). 

We measure the silicate strength Sio, 
„ , fobs(lO^im) 

SiQ = In , (1) 

/„,„,( 10/im) 

using the method of Spoon et al. (2007) where fobsiiOfim) 
is the observed flux density at the peak of the lO/zm feature 



and fconti^Ofxm) is the continuum flux at the peak wavelength, 
extrapolated from the continuum to either side. The silicate 
absorption depths are Sio= -1.0, -1.9, -3.6, and -2.2 for 
70Bootes5, 70Bootes6, 70Bootes8, and TOBootesll respec- 
tively. These correspond to very heavy absorption, similar to 
the median for ULIRGs in Hao et al. (2006) and Spoon et al. 
(2007) (Sio= -1.6) and places them among the most heavily 
absorbed sources known (the most extreme example known 
has Sio= -4.0; Spoon et al. 2006). Hao et al. (2006) find 
that the silicate strength correlates with the infrared slope at 
high mid-IR wavelengths (14-30/im) but not at lower mid-IR 
wavelengths (5-14 /im), supporting the idea that the silicate 
feature arises in the cooler dust (although note that no corre- 
lation between the f^(60/im)/f,y(25/xm) flux density ratio and 
silicate strength is seen by Imanishi et al. (2007) in their sam- 
ple of low redshift ULIRGs). Give that our selection criteria 
preferentially picks sources which are bright at 70 //,m, the 
large silicate strengths for the absorption-dominated sources 
are consistent with the findings of Hao et al. (2006). 

Are these sources powered by AGN or starburst activity? 
The individual spectra show no significant PAH-emission fea- 
tures (although there is a hint of the 7.7, 8.6, and 12.7/im PAH 
features in the averaged spectrum in Figure 2), so they appear 
most similar to class 3A (Uttle to no PAH emission, strong 
silicate absorption) in the scheme of Spoon et al. (2007). The 
IRS spectra are similar to that of 2 sources with deep silicate 
absorption strengths in a sample of 87 local Seyfert galaxies 
presented by Buchanan et al. (2006). Farrah et al. (2007) sug- 
gest that sources with very deep silicate absorption are likely 
to be AGN-dominated sources. The four sources have the 
largest 6/zm rest-frame luminosities in the sample and large 
infrared luminosities, again suggesting AGN activity. We ob- 
serve no significant [Ne V] emission line that we might expect 
for AGN-dominated sources. However, given that [Ne V] of- 
ten has a low equivalent width and that there are only 4 IRS 
spectra with low signal-to-noise at the observed wavelength 
of [Ne V] (see noise spectrum in Figure 1), this is not strong 
evidence against a large contribution from AGN emission. 

Because AGN-dominated galaxies typically have shallower 
mid-IR slopes than starburst-dominated galaxies (e.g., Brand 
et al. 2006), measuring flux density ratios at wavelengths 
free of strong absorption or emission lines may help in 
determining their primary power source. For heavily ab- 
sorbed sources, estimating the slope at large rest-frame mid- 
infrared wavelengths is important because dust may ab- 
sorb near-IR photons and re-emit them at longer wave- 
lengths, steepening the spectral slope. Figure 4 shows 
that the absorption-dominated sources have observed-frame 
i/f,y(70/im)/i^fi^(24/im) flux density ratios that are smaller 
than the PAH-dominated sources and closer to that of the 
X-ray selected AGN presented in Brand et al. (2007), sug- 
gesting that they have shallower infrared spectral slopes in- 
dicative of AGN-dominated sources. We note that one must 
be cautious in interpreting these results since the flux densi- 
ties are in the observed frame and the rest-frame flux den- 
sity ratio may vary with redshift. Although we can directly 
measure rest-frame flux density ratios at shorter wavelengths, 
these wavelengths are more likely to be affected by absorp- 
tion. The ft^(15/im)/f,y(6/im) flux density ratios of absorption- 
dominated sources are generally larger than that of the PAH- 
dominated sources. We suggest that strong dust absorption 
in the absorption-dominated sources results in a steeper in- 
frared slope at these shorter infrared wavelengths, despite the 
fact that they have small !/f!y(70/im)/z/f^(24/im) flux density 



ratios and are likely to be dominated by AGN. 

We conclude that these are examples of very obscured 
sources, which are hkely to be powered by AGN. However, 
we cannot rule out the possibility of them being powered by 
deeply embedded starbursts. If the mid-IR emission of these 
sources is dominated by AGN, this implies that ^36% of 
sources selected to have luminous 70 fim emission are AGN- 
dominated. Perhaps the 70 iim emission, which traces cooler 
dust, is also powered largely by AGN (see e.g., Imanishi et al. 
2007). In this cases, large 70 fim luminosities be generated 
because the AGN are so heavily absorbed by dust that a large 
fraction of the hot dust emission is absorbed and re-emitted 
at longer wavelengths. Alternatively, if a large fraction of 
the mid-IR emission originates from star-formation activity, 
this implies that star-formation activity as traced by mid- 
IR signatures can be hidden by obscuration and/or low level 
AGN activity. Follow-up multi-wavelength observations (op- 
tical spectroscopy, deep X-ray and radio imaging) are clearly 
needed to fully determine the nature of these rare sources. 

5. INFRARED AND INFRARED-TO-OPTICAL COLORS OF 70 fjM 

SELECTED SOURCES 

To help provide further insight as to the nature of the 
PAH-dominated and absorption-dominated sources, we show 
the IRAC color-color diagram for our IRS sources in Fig- 
ure 5a. For both the PAH-dominated and absorption- 
dominated sources, the colors are consistent with the expected 
positions of starburst-dominated relatively low redshift galax- 
ies (see also Yan et al. 2004; Sajina et al. 2005). Although 
we have concluded in previous sections that the absorption- 
dominated sources are probably largely powered by AGN 
emission, they do not fall into the "AGN wedge" of Stern et al. 
(2005). This is likely due to very heavy obscuration of the 
AGN. The IRAC colors of the absorption-dominated sources 
are consistent with their emission coming predominantly from 
heavily reddened star-formation processes (see reddening ar- 
row estimated from the Ri/=3. 1 dust model of Draine (2003)). 
The IRAC colors are redder than the PAH-dominated sources, 
consistent with them being more dusty and obscured than the 
PAH-dominated sources. 

We also plot the IRAC colors of the starburst-dominated 
sources selected by Weedman et al. (2006a) to have optical 
and near-lR SEDs showing a luminosity peak from stellar 
photospheric emission at 1.0 < z < 1.9. These sources have 
(by design) red [3.6]-[4.5] colors but blue [5.8]-[8.0] colors 
and confirm that sources in this region tend to be high red- 
shift starburst-dominated sources. This region becomes more 
heavily populated at lower 24 /xm fluxes, consistent with the 
population becoming more dominated by starbursts sources 
at fainter f24 as found by Brand et al. (2006). To inves- 
tigate the effects of obscuration, we show the 24-to-8 /xm 
vs. 24-to-0.7 /im color-color diagram for our sources in Fig- 
ure 5b, (see also Yan et al. (2004) for the location of differ- 
ent populations in this color-color space). The 24-to-0.7 /im 
color is a good indicator of obscuration whereas the 24-to- 
8 /im color can be used as a crude measure of the spectral 
slope and hence the dust temperature distribution and whether 
the source is AGN- or starburst-dominated (Brand et al. 2006 
find that AGN-dominated sources tend to have f(24:8)<0.3 
and starburst-dominated sources tend to have f(24:8)>0.3; 
see also Yan et al. 2004). The PAH-dominated 70 /zm se- 
lected sources have lower f(24:R) ratios than the absorption- 
dominated sources. Because PAH emission features can 
faU into the 8/xm IRAC wave-band for z <0.6, the 24-to-8 
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Fig. 5. — IRAC color-color diagram (top) and 24-to-8 fim vs. 24-to- 
0.7 /im color-color diagram (bottom) for 70 /im selected IRS sources from 
Bootes (large blue symbols) and starburst-dominated sources from Weed- 
man et al. (2006a) (smaller red squares). The Bootes sources are split into 
silicate absorption-dominated sources (large blue downwards triangles) and 
PAH-dominated sources (large blue squares). We show the distribution of all 
« 10,000 MIPS sources with (24 >0.5 miy in greyscale. The rough regions 
expected to be inhabited by different populations are also shown. The black 
dashed line in the top figure is that proposed by Stern et al. (2005) to em- 
pirically separate AGN from Galactic stars and normal galaxies. The arrow 
shows the reddening curve estimated from the Rv=3.1 dust model of Draine 
(2003) for a z=l source. In the bottom figure, the vertical dashed fine shows 
the R-[24] >14 (log(!/f„(24)M„(R))> 1.43) criteria used to selected the 
powerful obscured sources presented in Houck et al. (2005) and Weedman 
et al. (2006b) and the horizontal dashed line shows the 24-to-8/xm color cri- 
teria used by Brand et al. (2006) to roughly divide steeper spectrum starburst 
sources from shallow spectrum AGN-dominated sources at z >0.6. 

fim color is not a good indicator for most of these sources. 
The absorption-dominated 70/im sources have higher f(24:R) 
ratios than the PAH-dominated sources and steep spectral 
slopes, again consistent with them being more obscured. The 
starburst sources from Weedman et al. (2006a) have high 
f(24:R) ratios and f(24:8) ratios consistent with them being 
higher redshift, obscured starburst-dominated sources. 

6. MULTI- WAVELENGTH SEDS 

In Figure 6, we show the rest-frame multi-wavelength 
R, I, K, 3.6/im, 4.5/im, 5.8/im, 8/im, 24/im, 70/im, 160/im) 
SEDs for the PAH-dominated and absorption-dominated 
sources. We also show the SEDs of M 82, Arp 220, and Mrk 
231. M 82 is the nearest example of a star-bursting galaxy 
and is often considered the prototype of the starburst phe- 
nomena (e.g., Rieke et al. 1980). Arp 220 is a starburst- 
dominated ULIRG with strong dust extinction (e.g., Soifer 
et al. 1984; Sturm et al. 1996). Mrk 231 is a local ob- 
scured AGN with a luminous circumnuclear starburst (e.g.. 



Solomon et al. 1992; Weedman et al. 2005). The individ- 
ual SEDs are normalized to the same flux density at 3/im. 
All but one of the PAH-dominated sources are well fit by M 
82 blue- wards of 15/im, consistent with them being starburst- 
dominated sources. 70BooteslO (shown by the red points) is 
better fit by Mrk 231, suggesting that there is a large contribu- 
tion of AGN emission in this source. 70BooteslO is perhaps 
one of the most likely of the PAH-dominated sources to host 
an energetically important AGN, since it has the lowest 6.2/im 
PAH equivalent width of all the PAH-dominated sources. The 
far-IR flux densities are between that of M 82 and Arp 220, 
suggesting that the PAH-dominated sources have a large rel- 
ative cool to hot dust content that is between these two lo- 
cal sources. The absorption-dominated sources have a larger 
variation in their SED shapes. None of the sources exhibit 
a strong near-IR "bump" which is characteristic of starburst 
galaxies. The far-IR slopes are steeper than that of Mrk 231, 
suggesting that these sources have a larger cool to hot dust ra- 
tio than Mrk 23 1 . This may result from our selection criteria 
which requires the galaxies to be very luminous at 70/im. 

70um PAH sources 
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Fig. 6. — Optical to far-IR SEDs for PAH-dominated sources (top) and 
absorption-dominated sources (bottom). The individual spectra are scaled to 
have the same 3fj.m flux density. Also plotted are SED templates of Mrk231 
(solid puiple line), Aip220 (dashed blue line; Silva et al. 1998), and M82 
(dotted red line; Silva et al. 1998). For the PAH-dominated spectra, most 
SEDs are similar to M82 below ~15/^m. 70BooteslO is not weU fit by M82 
and its SED points are shown in red. 

7. CONCLUSIONS 

We have presented Spitzer IRS mid-infrared spectra of a 
small sample of 1 1 sources from the Spitzer MIPS survey of 
the Bootes field using selection criteria based on 70/im de- 
tections (fvo >30 mJy) and faint optical magnitudes (R > 20). 
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All the sources lie in the redshift range 0.3 < z < 1.3, im- 
plying very large mid- and far-infrared luminosities. The IRS 
spectra of the eleven galaxies show either PAH emission fea- 
tures (7/11) or deep silicate absorption (4/11) features. The 
IRS spectra of the seven PAH-dominated sources are typical 
of classical starbursts, but with extremely large mid-infrared 
luminosities. This sample contains sources with very high 
PAH luminosities, implying star-formation rates of up to 700 
Mq yr"^ The PAH-dominated sources tend to have lower 
z/f,^(I5/im)/i/fy(6/xm) flux density ratios than that of lower 
redshift starburst galaxies with lower infrared luminosities. 
A larger sample is needed to show whether this is consistent 
with a contribution to the f^(6/Ltm) luminosity from hot dust 
associated with an AGN in more luminous sources. 

The four absorption-dominated sources have very deep sil- 
icate absorption features that are larger than the median of 
that of local ULIRGs and their large !/f,y(15//,m)/i^f^(6/im) 
flux density ratios suggest that their mid-infrared spectral 
slope is steep. This suggests that the energy source is em- 
bedded in large volumes of cool dust. The very red [5.8]- 
[8.0] and relatively red |3.6]-[4.5] colors also suggest a large 
amount of dust along the line of sight to these sources. 
The deep silicate absorption, small j^f^(70/um)/i^fi,(24/im) 
flux density ratio (which is less affected by dust obscuration 
than estimates from lower infrared wavelengths), and multi- 
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wavelength SEDs suggests that they are likely dominated by 
AGN, although we cannot rule out that they are powered by 
very obscured starbursts. Follow-up multi-wavelength obser- 
vations such as optical spectroscopy, and deep X-ray and ra- 
dio observations may help to confirm their power source. 

We thank our coUeagues on tiie NDWFS, MIPS, IRS, 
AGES, and IRAC teams. KB is supported by the Giacconi 
fellowship at STScI. Support for this work by the IRS GTO 
team at Cornell University was provided by NASA through 
Contract Number 1257184 issued by JPL/Caltech. This re- 
search is partially supported by the National Optical Astron- 
omy Observatory which is operated by the Association of 
Universities for Research in Astronomy, Inc. (AURA) un- 
der a cooperative agreement with the National Science Foun- 
dation. This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology under a 
contract with NASA. The Spitzer /MIPS and IRAC surveys of 
the Bootes region were obtained using GTO time provided by 
the Spitzer Infrared Spectrograph Team (James Houck, P.I.), 
M. Rieke, and flie IRAC Team (G. Fazio, PL). IRAC is sup- 
ported in part through contract 960541 issued by JPL. The 
IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through 
the Jet Propulsion Laboratory and the Ames Research Center. 



Rieke, G. H., Lebofsky, M. J., Thompson, R. I., Low, F. J., & Tokunaga, A. T. 
1980, ApJ, 238, 24 

Rieke, G. H., Young, E. T, Engelbracht, C. W., Kelly, D. M., Low, E J., 

Haller, E. E., Beeman, J. W., Gordon, K. D., et al. 2004, ApJS, 154, 25 
Sajina, A., Lacy, M., & Scott, D. 2005, ApJ, 621, 256 
Sajina, A., Yan, L., Armus, L., Choi, P., Fadda, D., Helou, G., & Spoon, H. 

2007, ArXiv e-prints, 704 
Sanders, D. B., Soifer, B. T, Ellas, J. H., Madore, B. E, Matthews, K., 

Neugebauer, G., & Scoville, N. Z. 1988, ApJ, 325, 74 
Saunders, W., Rowan-Robinson, M., Lawrence, A., Efstathiou, G., Kaiser, 

N., Elhs, R. S., & Frenk, C. S. 1990, MNRAS, 242, 318 
Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 103 
Smith, J. D. T, Draine, B. T., Dale, D. A., Moustakas, J., Kennicutt, Jr., R. C, 

Helou, G., Armus, L., Roussel, H., et al. 2007, ApJ, 656, 770 
Soifer, B. T., Neugebauer, G., Helou, G., Lonsdale, C. J., Hacking, P., Rice, 

W, Houck, J. R., Low, E J., & Rowan-Robinson, M. 1984, ApJ, 283, LI 
Soifer, B. T., Sanders, D. B., Madore, B. F., Neugebauer, G., Danielson, G. E., 

EUas, J. H., Lonsdale, C. J., & Rice, W. L. 1987, ApJ, 320, 238 
Solomon, R M., Downes, D., & Radford, S. J. E. 1992, ApJ, 387, L55 
Spoon, H. W. W, Marshall, J. A., Houck, J. R., EUtzur, M., Hao, L., Armus, 

L., Brandl, B. R., & Charmandaris, V. 2007, ApJ, 654, L49 
Spoon, H. W. W, Tielens, A. G. G. M., Armus, L., Sloan, G. C, Sargent, B., 

Cami, J., Charmandaris, V., Houck, J. R., & Soifer, B. T. 2006, ApJ, 638, 

759 

Stem, D., Eisenhardt, P, Gorjian, V., Kochanek, C. S., Caldwell, N., 

Eisenstein, D., Brodwin, M., Brown, M. J. 1., et al. 2005, ApJ, 631, 163 
Sturm, E., Lutz, D., Genzel, R., Sternberg, A., Egami, E., Kunze, D., 

Rigopoulou, D., Bauer, 0. H., et al. 1996, A&A, 315, L133 
Tran, Q. D., Lutz, D., Genzel, R., Rigopoulou, D., Spoon, H. W. W, Sturm, 

E., Gerin, M., Hines, D. C, et al. 2001, ApJ, 552, 527 
Weedman, D. W., Hao, L., Higdon, S. J. U., Devost, D., Wu, Y, 

Charmandaris, V., Brandl, B., Bass, E., & Houck, J. R. 2005, ApJ, 633, 

706 

Weedman, D., Polletta, M., Lonsdale, C. J., Wilkes, B. J., Siana, B., Houck, 

J. R., Surace, J., Shupe, D., et al. 2006a, ApJ, 653, 101 
Weedman, D. W., Soifer, B. T, Hao, L., Higdon, J. L., Higdon, S. J. U., 

Houck, J. R., Le Floc'h, E., Brown, M. J. I., et al. 20()6b, ApJ, 651, 101 
Yan, L., Chary, R., Armus, L., Teplitz, H., Helou, G., Frayer, D., Fadda, D., 

Surace, J., & Choi, P 2005, ApJ, 628, 604 
Yan, L., Helou, G., Fadda, D., Marleau, R R., Lacy, M., Wilson, G., Soifer, 

B. T, Drozdovsky, I., et al. 2004, ApJS, 154, 60 
Yan, L., Sajina, A., Padda, D., Choi, P., Armus, L., Helou, G., Teplitz, H., 

Frayer, D., & Surace, J. 2007, ArXiy Astrophysics e-prints 



